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Abstract Molecules in comets are formed through chemical oxidation 
reactions induced by radiation. Thesereactions can be simulated in laboratory 
experiments applying gamma radiation to samples at low temperatures. The 
kinetics of the induced reactions can be modeled by a system of coupled non¬ 
linear ordinary differential equations describing the mass balance of all of 
the species involved. However, finding a traditional solution to this system 
is difficult because of the large number of reactions involved, the need to 
solve all of the equations simultaneously, and the strong dependence on the 
initial conditions due to the non-linear character of the equations. For each 
species, the mass-balance equation includes all of the reaction rates leading to 
production (source terms) and to destruction (sink terms). In this sense, each 
equation is analogous to the prey-predator model, with the sink terms consider 
to be the “prey” and the source terms as the “predators”. Due to this, we can 
use an agent-based model to follow the kinetics of the chemical reactions. 
In this paper, we present a code in Python for an agent-based model of the 
chemical oxidation of ferrous ions (Fe 2+ ) induced by gamma radiation and 
in the presence of molecular oxygen. We compare the results that this code 
produces for molar concentrations of Fe 3+ over time with those obtained in the 
laboratory. 

Keywords: Chemical reactions; radiation; kinetics of reactions; prey-predator 
model; Fricke-model. 
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1. INTRODUCTION 

Chemical oxidation reactions induced by radiation at low temperatures 
are fundamental in laboratory experiments that simulate the formation of 
biologically relevant molecules in icy bodies l ik e comets [4,7, 10, 11]. In these 
simulated experiments, dose measurements can indicate the amount of energy 
deposited by the gamma radiation on the samples providing a quantitative 
measure of the changes induced by ionizing radiation. Chemical analysis, 
pulse radiolysis, and competition kinetic methods have been employed 
experimentally to obtain the large number of rate constants involved in the 
reactions as in [3, 6]. 

To describe the products generated by the interactions among the different 
molecules and their chemical reactions under radiation, mathematical models 
exist that describe the kinetics of the chemical reactions as coupled nonlinear 
ordinary differential equations (NODEs) [5, 15, 16]. Solving these models 
is difficult because ofthe large number of reactions involved, the need to 
solve all of the equations simultaneously, and the strong dependence on 
initial conditions due to the non-linear character of the equations [5,9]. To 
overcome these numerical troubles, it has been used Monte-Carlo simulations 
[17]. Here, we propose an agent-based model to simulate the chemical 
evolution of oxidation reactions among ferrous ions under radiation. This is 
a modified version of the prey-predator model [2, 8, 13, 14], in which each 
chemical species involved is considered to be an agent that can interact with 
other species with known reaction rates. The agent-model is programmed 
in Python. The radiation is considered as a factor that promotes product 
creation/destruction. This original model can reproduce the formation of 
products and the consumption of ferrous and oxygen ions in a laboratory 
reaction [15]. 


2. EXPERIMENTAL SETUP FOR THE OXIDATION OF 
FERROUS IONS 

We irradiated oxygenated solutions of iron salt (1.6*10" 3 M Fe,S0 4 -7H,0) in 
sulfuric acid 0.4M in a Cobalt 60 source at a dose rate of 0.93Gy/min (following 
the setup of [10]). The samples were in Pyrex cells surrounded by a Dewar 
flask filled with different cooling agents. In the frozen case, the samples were 
melted after irradiation, and the temperature was adjusted to 25°C prior to the 
spectroscopic measurement of the Fe 3+ concentrations. Our first approximation 
was with the data at 295 K. 
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3. OXIDATION REACTION MODEL 

3.1. Chemical oxidation reactions: Ferrous ions(Fricke dosimeter) 

In the case of the Fricke dosimeter, the chemical oxidation reactions for the 
Fe 2+ are well known [15]: 
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H°+0 2 ^HO° 

(i) 

Q ~aq + H 3 0 + —> H° 

(2) 

Fe 2+ + OH Fe 3+ + OH" 

(3) 

Fe 2+ + H 2 0 2 ^ Fe 3+ + OH + OH“ 

(4) 

Fe 2+ + H0 2 ^ Fe 3+ + H0 2 

(5) 

HO"+H + %H 2 0 2 

(6) 

h 2 o 2 4 ho;+h + 

(7) 


In this case, there are seven chemical reactions that involve 11 species (Fe 2+ , 
Fe 3+ , H 2 0, H", HO,', HO,, 0 2 , H 3 0 + , OH, OH - , and H 2 0 2 ), and seven known 
reaction rates (r). We considered H 3 0 + , Fe 2+ , O,, and H,Oto be compounding 
the pre-irradiated solution, so initially there are no other species (the initial 
concentrations of Fe 3+ , H°, HO,’, HO,,', OH, OH', and H,0, are zero). We used 
the standard concentrations for the Fricke dosimeter (from 0.1 to 0.001M for 
Fe,S0 4 -7H 2 0), in which the oxidation of ferrous ions takes place in an acid 
medium (H,S0 4 , 0.4M). The concentration of oxygen is based on its solubility 
in water, as given by [12]. Then, the values for the initial conditions were 0.0016 
M Fe 2+ , 0.0120 M O^, and 0.4000 M H 3 0 + .The reaction rates were:r 1 ~2.0*10 10 , 
r~ 2.3* 10 10 , r 3 ~3.0*10 8 , r 4 ~50, r~ r~ 7.3* 10 5 , and r 7 ~10 5 in M 's 1 . 
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3.2. Mathematical model 

We are interested in monitoring the concentration of the 11 species involved 
as the experiment progresses, to follow the kinetics of the chemical reactions 
(Eqs. 1-7) induced by the radiation. These chemical reactions can be modeled 
through differential equations that establish the mass balance for each species 
involved. The mass balance equations predict the expected change in the 
molar concentration of the species over time, and take into account all of the 
reaction rates leading to production (source terms) and to destruction (sink 
terms). In this sense, it is analogous to the prey-predator model, withthe sink 
terms being considered the “prey” and the source terms as the “predators”. 
Importantly, a species can be a prey in one of the mass balance equations; 
and in another equation, the same species can play the role of predator. In 
this model, we consider that each balance equation reflects the interaction of 
at most 2 species, and that such a reaction, whatever its rate, must lead to the 
production or consumption of a certain species. In general, the NODEs are of 
the form given by [5]: 

< 8 ) 

ul j k 

m 

where ^ (/) denotes the molar concentration of the / -species at time t, ft is 
the number of molecules of the i -species produced i-n the radiolysis process 
induced by the gamma radiation, and r‘\ denotes the reaction rate constant for 
the two reactive species j ,k that produce the specie (z) .The radiolysis term, 
fi, is an external source in the balance equations (8), which we propose as a 
linear function of the radiation dose intensity I d , given in rad/h, and of the 
temperature T, as in [5]: 




6.2*10" M. 


3.61V. M 


G,\(T)l d , 


H.O 


where N A is the Avogadro number (6.022xl0 23 molecules), M i is the 
molecular mass of the species/, M H O is the molecular mass of water (18.02g/ 
mol), and G i is the radiochemical yield that gives the amount of species/ 
produced when lOOeV are absorbed by the system, in this case 15.6 molecules 
of Fe 3+ . 
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The system of chemical reactions given by Eqs. (l)-(7) can be rewritten as 
a system of 1 INODEs, Eqs. (8), one for each species involved.For the ferrous 
ion, Fe 2+ , the NODE is: 
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(0 

dt 


~ / Fe 2+ {l<i’T) ~ r 3 N OH (i)iV Fe 2 + (t) 

-'4^H 2 0 2 (O^V (0 “^HO; ( *) N Fe 2 + (0 > 


(9) 


where only Fe 2+ is consumed in agreement with Eqs. (3), (4) and (5). This Fe 2+ 
is responsible for the production of Fe 3+ , through equation: 


dNf dt^ =f ^ (^’ r ) +r ^oH(0^V (0 («>) 

+ r 4^H 2 0 2 (O^V (0 + r 5^ HO ° (O^V (0 ' 

If we consider the radical hydrogen, H", as the species, the corresponding 
NODE is: 

dr ^r=f*- ( I - T ) +r i N ,„ (')V* ('HN (ON ('). CD 

where the production terms involve Eq. (2), and the consumption terms come 
from eq. (1). Notice that H°is produced at the same rateas H 3 0 + and e~ are 
consumed. In the case of oxygen, 0 2 , the corresponding equation only has a 
loss term that arises from eq. (1): 

dN 0 (t ) . . . . . . 

—| ii = /o.(4.r)-';N,(0NW- <> 2 > 

If we continue doing this, we get a system of 11 coupled NODEs, which can 
be solved numerically through traditional methods, as the one's presented in 
[16], or by considering an agent-based model [1], 


4. NUMERICAL SOLUTION 

A numerical solution for the NODEs Eqs. (8) for all of the species can be 
found by transforming each differential equation to a difference one: 
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j k 


+N i (t- 1 ) 


1 


(13) 


For example, the molar concentration of Fe 3+ at moment t can be obtained from 
the molar concentrations of OH, Fe 2+ , H,0„ and 110 at the previous time t- 1 
through: 


A... W="*»('- 0 +/„» M+r,N m (i- 1 ) ((- 1 ) 

+'^ h i0 , ('-O + '-sV (<-l)iV F ,,.('-l) ■ (14) 


and Eq. (11) for H takes the form: 

The coupled system of eqs. (13) can be solved using an agent-based model 
with sources and sinks that can be interpreted as prey and predators. We made 
a Python program to follow the molar concentrations of each species. The code 
initializes a region of 1,000* 1,000 pixels, with 4,000 randomly allocated points 
associated with H 3 0 + , 120 points of 0 2 , and 16 of Fe 2+ (see Figure la). Each 
position is checked during the next time instant, if two neighbor species can 
react in accordance to eqs. (1) to (7), the reaction takes place consuming the 
“prey” elements by the “predators” produced. The rate of reaction is determined 
by the rj\ . After actualizing the entire net, the procedure is repeated until the 
steady-state is reached (see Figure lb). In the program we can follow the 
molar concentration of each species as a function of time. In Figure 2, the 
molar concentration of Fe 3+ is plotted and compared with experimental data. 

One of the common problems in following the kinetics of many chemical 
reactions is that the reaction rates are unknown, especially at low temperatures. 
Our agent-based model can be used to evaluate some of these rates using 
experimentally determined concentration values. For this, we ran the model 
with different values of unknown rates until we reached the best experimental 
approximation. In Figure 3, we plotted the relative error in the evaluation of 
the different rates when we did not know from one to four rates. Our model 
had an evaluation error of less than 10% when we did not know the values of 
one or two parameters. The worst approximation was given to r 4 , the smallest 
reaction rate. 
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Figure 1: Agent-based spatial distribution of the species (a) initially and (b) at steady- 
state. 



time 


Figure 2: Molar concentration of Fe 3+ . The points correspond to experimental values 
and the continuous line to the agent-based results, at room temperature. 
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Figure 3: Evaluation error of the reaction rates obtained from the agent-based model. 
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5. CONCLUSION 

In this paper we described the kinetics of the oxidation of Ferrous ions induced 
by gamma radiation at room temperature in a Fricke dosimeter. Instead of 
solving the system of coupled non-linear ordinary differential equations 
that describe the mass balance of all the species involved (Eqs. 8), here we 
implemented an original agent based model that includes in Eqs. (13) all of 
the reaction rates as source and sink terms (like in a prey-predator model). 
The agreement of the results from our agent-based model developed in Python 
with the experimental results suggests that the numerical model is reliable for 
analyzing the kinetics of chemical oxidation reactions induced by radiation. 
This model is very simple and robust, and maybe could be applied to different 
chemical reactions. We expectthat it can be used to understand the kinetics of 
reactions induced by radiation when the temperature of the irradiated medium 
decreases below the freezing point. Further studies and experiments are needed 
to reveal the utility of our model. 
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